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Abstract 

The survival probability of immobile targets, annihilated by a population of random walkers on 
inhomogeneous discrete structures, such as disordered solids, glasses, fractals, polymer networks 
and gels, is analytically investigated. It is shown that, while it cannot in general be related to the 
number of distinct visited points, as in the case of homogeneous lattices, in the case of bounded 
coordination numbers its asymptotic behaviour at large times can still be expressed in terms of 
the spectral dimension d, and its exact analytical expression is given. The results show that 
the asymptotic survival probability is site independent on recurrent structures (d < 2), while on 
transient structures {d > 2) it can strongly depend on the target position, and such a dependence 
is explicitly calculated. 

PACS numbers: 82.33.-z, 82.39.Rt, 05.40.-a 
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The kinetics of diffusion limited reactions is deeply affected by geometry and topology, 
and their description in terms of random walks models is a very powerful tool to explore 
such a dependence In some cases, especially on regular lattices, it is possible to 

establish simple analytical relations between survival probabilities of chemical species and 
basic random walks functions, making explicit the dependence on universal geometrical 
parameters such as spatial dimension. The case of irregular structures is more complex, due 
to the absence of symmetries allowing for a reduction of the degrees of freedom involved 
in analytical calculations. Therefore, in general, it is not possible to relate time decays 
and survival probabilities to simple geometrical parameters characterizing the underlying 
structures. 

In this paper, we face a well known kind of reaction, namely the A+B —>■ B process, where 
the A species is immobile, which is known as the target reaction (A is called the target) j^. 
This target reaction model is usually introduced to describe a variety of physical, chemical 
and biochemical processes, such as the Williams- Watts dielectric relaxation in polymers and 

n n 

glasses as well as the poisoning of surface catalysts and of immobilized enzymes [5[. 
Such a target decay problem has been extensively studied analytically, since it admits exact 
solutions on homogeneous structures, i.e. on structures where all sites are topologically 
equivalent Js, 6]. In these cases, it can be shown that the target survival probability decays 
as a negative exponential of S(t), the mean number of distinct sites visited by a random 
walk after time t. The case of inhomogeneous structures has been recently investigated 

Q 

for some particular networks, namely Small World Networks (SWN) [71] and Scale Free 
Networks (SFN) 8|], showing that this simple dependence on S(t) no longer holds, and that 
the deviations from such a behavior are particularly evident at large times. These networks, 
besides being inhomogeneous, exhibit a peculiar feature: the coordination numbers are not 
bounded from above in the thermodynamic limit: in other words the maximum coordination 
number diverges for N ^ oo, N being the number of sites. 

Many interesting real inhomogeneous system, however, don't present such a property and 
the maximum coordination number is finite even for — > oo : this is the case, for example, 
of disordered solids, glasses, fractals, polymer networks and gels. All these structures can be 
mathematically described in terms of "physical graphs" , and the applications of ideas and 
techniques of algebraic graph theory allows us to obtain analytical results in spite of the lack 
of invariance and symmetry jol] . In the following we present the mathematical formulation of 
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the target decay problem on physical graphs, and we obtain an exact analytical expression 
for the target survival probability at large times, showing that it can be expressed as a 
negative exponential of t for transient graphs, and of for recurrent graphs, d being the 
graph spectral dimension. ^ 
Let us begin with some basic definitions [9j]: 

A graph ^ is a countable set V of vertices (or sites) (i) connected pairwise by a set E of 
unoriented links (or bonds) = (j, 0- 

The graph topology can be algebraically represented introducing its adjacency matrix Aij 
given by: 



and the coordination number of site i, which is the number of nearest neighbors of i, is given 

by Zi = T^j^ij- 

The discrete time simple random walk on a graph Q is defined by assuming that at 
each time step t the walker can only jump to a nearest neighbor site, and that all nearest 
neighbor sites can be reached with the same probability. Therefore, we can define the 
jumping probabilities pij between sites i and j by 



where Zij = Zi6ij. 

Now we introduce the functions Pij{t), each representing the probability of being in site 
j at time t for a walker starting from site i at time 0, and the first passage probabilities 
Fij{t), each representing, for j ^ i, the conditional probability for a walker starting from i 
of reaching for the first time the site j in t steps, and, for i = j , the probability of returning 
to the starting point i for the first time after t steps (-Fm(O) = 0). 

The fundamental relation between the Pij{t) and the Fij{t) is given by 

t 

P^At) = E P^A'')Pnit -r)+ S.M (3) 

r=0 

Introducing the generating functions -Pij (A) and Fij{X) by the definition 




(1) 



Pij 




(2) 



oo 




(4) 
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from eq. ([3]) we obtain the simpler relation 



4(A)=F,,(A)P,,(A) + 5., (5) 
which will be useful in the proof of our main results. 

On infinite graphs, representing real systems in the thermodynamic limit, Pii{t) vanishes 
for t ^ oo. If the graph can be embedded in a finite dimensional Euclidean space, and 
if the coordination numbers are bounded, i.e. if 3 Zmax I < z^ax^i ^ V, then Puit) 
vanishes tipically as a power law, whose exponent allows to define the so called (local) 
spectral dimension which is the natural generalization of the Euclidean dimension for 



dynamical processes 
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Puit) r^Pit-"^^^ for t-^oo \fieg (6) 



with d > 1, where Pi = po Zi (i.e. it depends only on the coordination number of i) if 
d < 2. For d < 2, Fij{l) = 1, i.e., the probability of ever reaching any site starting from any 
site is 1, and the graph is called recurrent. For d > 2, -Fij(l) < l,and the graph is called 
transient. 

Now we can define the target decay problem on an infinite graph, following the formalism 
introduced in [3]. 

Since each target decays independently of the other ones, we can study the decay of a 
single target without loss of generality. Let us suppose a target molecule A is placed at site 
k, while, at time t = 0, the B molecules are randomly and independently distributed over 
the other sites, with average site occupation number q. The occupation number distribution 
at each sites turns out to be Poissonian, and the probability p{n) of finding exactly n B 
molecules at a given site is p{n) = 2-^. For t > 0, the B molecules are moving randomly 
and independently according to the jumping probabilities ([2]), and the target A is annihilated 
when it is reached by one of them. Under these hypotheses, it has been shown that the 
survival probability ^k(t) of target A at time t is given by 



$;,(t) = e-''®''^ (7) 

where 
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t 



(8) 



Now, on homogeneous graphs, Fjfc(t) = F^iit), and, due to this symmetry, Qkif) = 
S{t) — 1, where S{t) is the number of distinct sites visited by a walker after t steps (on 
homogeneous graphs it is independent of the starting site k). This result gives rise to 
the well known results obtained on d-dimensional Euclidean lattices, where, for t — > oo, 
Qk{t) ~ for c/ = 1, Qk{t) ~ t/\nt for d = 2, and <S)k{t) ~ t for c/ > 3 Q. 

On inhomogeneous graphs, Fik{t) ^ Fki{t), and the simple relation mentioned above no 
longer holds, giving rise to a more complex behavior j?]. 

Let us proceed to the calculation of the asymptotic behavior of Qk(t) by introducing its 
generating function: 

oo ^ 

e,(A) = ^ A*e,(t) = ^^F,,(A) (9) 



t=0 i^k 



From ioY i ^ k, we get 



therefore 



MX) = = (10) 



Zk Ph{X) ^ p^^^^^ ^ ^Ph{X) ^^^^ 



Zmax Pkk{X) ^mmPfeA:(A) 

where z^^^ > 1 is the minimum coordination number. 

Moreover, since Pkiif) = 1 for every t, we have that ^ Pki{X) = (1 - A)~^ - Pkk{X) 
Therefore 



^ ^ <e.(A)<-^( ^^ -—^] (12) 



L - A)2Pfc,,(A) 1 - a; - ' ' - z^^ V(i - xyPkkix) 

Now we can proceed to the singularity analysis of QAX) in order to obtain the asymptotic 
behavior of Ofc(t) by applying Tauberian theorems 
From ([6]), we have 



PkkW 



PoZkT (^l-lj(l-A) 
Po Zk log(l - A)"^ 
Pkk (1) 



for < 2 
for d = 2 
for d> 2 



(13) 



therefore, from ([T^ . 



1 



1 



ZfcPO (l-A)2log(l-A)-i 
Zk 1 



for < 2 

for d = 2 
for d > 2 



(14) 



where ^ 



{ ZkPkk (1) (1-^)' 

Um (1 - \)Y.PkiW-t- = 1™ Y.PkiWj-/Y.PkiW^ with Z^in <Zk<Z„ 



is the weighted average of the inverse coordination numbers 1/zi, with weights Pki{X), for 
A ^ 1^. Notice that the hypothesis of boundedness of Zi could be replaced with the weaker 
condition Zk < oo, leaving our results unchanged. On recurrent graphs, since lim -^4tt = 1; 
uniformly in i, Zk turns out to be site independent: Zk = z. The singularities in A = 1 in 
( HM finally give us the following asymptotic behaviors: 



zpQ 7r| 



■t2 



t—*oo 



^ for 

zpo log t 



Zk 



, ZkPkk (1) 



for d < 2 
2 

t for d> 2 



(15) 



These results deserve some comments. First of all notice that, if the graph is homo- 
geneous, we exactly recover the usual asymptotic form of S{t): this happens not only for 
Euclidean lattices jsj, but also for regular ultrametric spaces {g], and for Bethe lattices, for 
which the spectral dimension turns out to be infinite pj|. In other cases, such as fractal 
structures, even if 9fc(t) is different from Skit), it turns out to have the same t dependence 
for t oo. 

For completeness sake, we notice that, in some particular cases, it is possible to have 
logarithmic corrections to the asymptotic behavior described by [12]: these corrections 
give rise to logarithmic corrections in ( ITSl) too, which are rather simple to calculate; we have 
neglected them only for simplicity's sake. 
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Moreover, for recurrent graphs, the asymptotic behavior of Qk{t) is site independent 
(0fc(t) even if the structure is inhomogeneous and the sites are not equivalent. 

This means that, in spite of inhomogeneity, the recurrent nature of the structure gives rise 
to asymptotic survival probabilities which are independent of the target position: this is the 



case, e.g. of many deterministic fractals such as the Sierpinski gasket and the T-fractal 
and also of recurrent random graphs generated with the constraint of bounded coordination 
number [l3 |. 

On the other hand, when the graph is transient, Qk{t) and are site-dependent even 

for t oo. Such a dependence, that only concerns the coefficient of the power of t, while 
the exponent is the same for all sites, is rather intriguing. Indeed, it is a peculiar feature of 
inhomogeneous transient graphs. In fact, only on these structures the probability fi of ever 
reaching site i, averaged over all possible starting sites, can depend on i. In other words, 
there can exist sites which are more likely to be visited than others, even at large times, and 
such a property gives rise to the asymptotic site dependence of 

From the point of view of the applications, the site-dependence is quite interesting, since 
it means that the target position can affect its survival probabilities at large times, and, in 
some cases, it is possible to know which sites are more likely to survive (notice that P^k (1) 
is the average number of visits to k for a walker starting from k itself). It is particularly 



noteworthy the case of the so called recurrent on the average transient graphs 



14| . such as. 



e.g., NTe) ("Nice Trees of dimension D") 151], where Pkk (1) is always finite but unbounded 



from above, i.e., for every r G R, it exists some k such that Pkk (1) > f^- In the case of the 
NT^), which are trees with branches of unbounded length, Pkk (1) is greater for points lying 
in longer branches: therefore a target placed in such sites has a greater asymptotic survival 
probability. 

The detailed investigation of these aspects is fundamental to design optimal reaction 
strategies based on geometry, as well as to understand target decay processes in complex 



biological systems [16 1 
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